Field-induced disorder in a gapped spin system with non-magnetic impurities 
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We study S = | dimers which are weakly coupled by three-dimensional antiferromagnetic inter- 
actions J'. It is commonly known that doping with non-magnetic impurities immediately induces 
a long-range Neel order. We show that application of an external magnetic field H may drive the 
system back to a disordered phase. We discuss the zero temperature phase diagram in the (H, J') 
plane and we propose suggestions for experiments. 

PACS numbers: 74.20Hi, 75.10Lp, 71.10+x 



r/3 



I 

C 

O 
o 



> 

m 
in 
m 
O 

o 

03 



o 
o 



X 



In the last decade, much attention has been paid to the 
physics of gapped spin systems doped with non-magnetic 



impurities. Typical examples include S = 4 ladders ElE 




J3j £ = | dimerized chains |6|, Haldane spin chains 
various two-dimensionally coupled systems H 0, 
llll I12| and systems of S = | dimers weakly coupled 
by three-dimensional (3D) interaction [l3L It was 

shown P, 0, @ that the presence of impurities induces 
the formation of a gapless continuum of low-lying states 
within the gap, which leads to long-range Neel ordering 
in the presence of arbitrarily small 3D interactions. 

On the other hand, considerable interest has been de- 
voted to the behavior of gapped spin systems in strong 
magnetic fields. When the field is strong enough to close 
the spectral gap, the system enters a new phase. In ab- 
sence of anisotropy, this new phase is critical in the purely 
one-dimensional (ID) case and exhibits long-range order 
in the presence of a weak 3D coupling. Experimentally, 
field-induced ordering has been studied for several sub- 
stances, one of the better understood examples being the 
dimer material TlCuCl 3 E1E1E3E3 

Consider an interacting 3D system of spin-^ dimers 
described by the following Hamiltonian: 

H = j S r> i ■ S ri 2 — h S^, a 

r ra 

+ J' J]] (Sr,i ' S>',2 + Sv,2 • Sy'.i)) (1) 

{rr') 

where J > and J' > are the intra- and interdimer 
exchange couplings, respectively, and H is the external 
magnetic field directed along the z axis. The vector r 
labels the dimers located at the sites of a lattice with 
coordination number Z, and (..) denotes summation only 
over neighboring dimers. In Q we have assumed that 
the couplings are unfrustrated so that each spin S r ^ a can 
be classified as belonging to one of the two sublattices 
(respectively, a = 1 or a = 2). 

In the absence of impurities and external field and for 
sufficiently weak interdimer coupling J' -C J the system 
has a singlet ground state without any magnetic order, 
and a finite gap A ~ J to the lowest excitation which 



is a triplet; known examples of materials exhibiting this 
type of ground state are KCuCl 3 EH and TlCuCl 3 HJ. 
If the 3D coupling J' exceeds some critical value (typi- 
cally of the order of J/Z), long-range Neel order appears; 
this situation is realized e.g. in NH4CUCI3 j2l|. An ex- 
ternal magnetic field closes the gap at the critical field 
H = H c = A, inducing a phase transition. At H > H c 
a finite magnetization along the field direction appears, 
accompanied by a staggered order in the plane perpen- 
dicular to the field; the U(l) symmetry in the xy plane 
is spontaneously broken. This transition may be viewed 
as the Bose- Einstein condensation of magnons EH • 

The presence of non-magnetic impurities drastically 
changes the above picture Q, 0, E3 • Impurities gener- 
ate unpaired spins, which develop an effective interaction 
with each other mediated by the intact dimers between 
them. According to the Lieb-Schulz-Mattis theorem, the 
interaction between two unpaired spins is antiferromag- 
netic if they belong to the same sublattice and ferromag- 
netic otherwise, so in total the unpaired spins encourage 
an antiferromagnetic ordering which at T = occurs at 
any arbitrarily small concentration of impurities c. Un- 
paired spins thus form a separate subsystem with the 
energy scale Jbr its dynamics set by the average effective 
interaction J -C A. The excitation spectrum of the im- 
pure system is determined by the gapless continuum of 
the low-lying states resulting from the ordered unpaired- 
spin subsystem; this continuum coexists with the states 
above the spin gap A of the pure system which survive 
with reduced weight. 

So, taken separately, both doping with non-magnetic 
impurities and application of an external field tend to 
induce ordering in a gapped spin system. The aim of 
the present Letter is to show that application of an ex- 
ternal field to a doped system may drive it back into 
the disordered phase, causing reentrant behavior of the 
Neel order as a function of H . The qualitative form of 
the T — phase diagram in the (H, J') space will be 
obtained from a mean-field treatment of the interaction 
between the dimer and unpaired spin subsystems. 

Dynamics of an intact dimer system.- To derive the 
effective interaction between two unpaired spins, we have 
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to consider the dynamics of the pure dimer system first. 
For this purpose we use the dimer field theory which 
may be viewed as a continuum version of the bond boson 
operator approach |l7| . The quantum state is formulated 
as the product of dimer states of the form 

H>) = (1-A 2 - B 2 ) l ' 2 \s) + + iB^tj), (2) 

j 

where \s) denotes the singlet state and \tj), j — x,y,z are 
the three triplet states of the spin dimer in the Cartesian 
basis. The vectors A, B are connected with the mag- 
netization M = (Si + S2) and sublattice magnetization 
L = (Si — S2) of the spin dimer as follows: 

M = 2(AxB), L = 2(1 - A 2 - B 2 ) 1/2 A . (3) 

The Lagrangian density (per dimer) can be written as 

C = 2hB ■ ^ - E dim - ^{(Vi) 2 + (VM) 2 } 
7 V 

E dim = J(A 2 + B 2 ) H — -j— (M 2 - L 2 ) — HM Z , (4) 

where we use the notation (VA) 2 = ^^dA/dxi) 2 . As- 
suming A, B <C 1, at the quadratic level one easily 
obtains the dispersion of a magnon with S z = p, as 
e(k) — A(l + k 2 ^ 2 ) 1 / 2 - pH, where 

A=[J(J-ZJ^, -(^) V2 - (5) 

Here £ is the spin correlation length, £ <C 1 if the 3D 
interaction is small and Z J' <C J (the lattice constant is 
set to unity). 

Interactions between unpaired spins. - The effective in- 
teraction between two unpaired spins can be estimated 
in second-order perturbation theory, 

J e(k) rvt, \£/ 

where K, is an unknown constant of the order of unity, de- 
termined by the matrix elements of the perturbation, r is 
the vector connecting the two spins, Ki is the MacDon- 
ald function, and we omit the oscillating sign depending 
on whether or not the spins belong to the same sublat- 
tice; we have used the H = expression for the magnon 
energy since, as will be clear below, we are interested in 
fields which are small comparing to the gap A. If r = 1, 
there is an additional direct (first-order in the perturba- 
tion J') interaction between the spins. Since ^< 1, one 
can write down the effective interaction as 

Jeff(r) = J 5 rA + Jir- 3 ^ 2 e- r ^, (7) 

where J - J' and J x = (2^/2i:f/ 2 KJ l (Z 5 J l /J) 1 / i . 
Since J e g is decaying very fast with increasing r, it is 
sufficient to take into account interactions only between 
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FIG. 1: A schematic mean-field T — phase diagram of 
the model Q.The dashed line corresponds to a pure system 
and the solid line to a finite concentration of impurities. The 
boundaries of the disordered and saturated phases are deter- 
mined by 11311 and 1141 1. respectively. 

those unpaired spins which are nearest neighbors. For a 
given impurity concentration c, the distribution of dis- 
tances R from an impurity to its nearest neighbor can be 
estimated in the continuum approximation (i.e. neglect- 
ing the presence of the lattice) as 

p(R) = c exp{ — — i? 3 }, / p (R)4ttR 2 dR = 1. (8) 

The distribution of effective exchange couplings J be- 
tween the neighboring unpaired spins is given by 

P(J) = J2p( R ) 5 (J- J es(R))- 
R 

Replacing the sum by an integral and substituting the 
exact value of p(l) = Zc, one obtains 

P{J) = Zc5(J - J ) + Preg(J), (9) 

P rcg (J) = 0( Jie~ 1/? - J) R\ exp{-4TCi??/3}, 

where Ri ~ £ ln(Ji/J) and is the Heaviside func- 
tion; note that P(J) has an upper cutoff. The peak 
of the regular part P rcg occurs at J = J poa k — 
Ji exp{— (47rc£ 3 ) -1 / 2 }, and the "regular average" is 

(J)reg = / JPreg(J) dJ ~ 4^Jie~ 1/ « > J pcak . (10) 

The interacting system of unpaired spins and intact 
dimers.- We go on to study the static properties (e.g. 
the sublattice magnetization) of the impure system in a 
magnetic field in a mean-field approximation by minimiz- 
ing the full energy E = E d im + E unp + E int . Here, -Edim 
is the energy of the intact dimers, Eq. (@J, E unp is the 
energy of the subsystem of unpaired spins and E{ nt the 
interaction energy between the unpaired spins and the 
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intact dimer subsystem. The unpaired spin subsystem is 
ordered, and we therefore describe it as a two-sublattice 
antifcrromagnet. Denoting the angle between the sublat- 
tice magnetization and the z axis as 9 and the average 
spin length of an unpaired spin as s , we write the inter- 
action energy (per dimer) as 

E int = cZJ's(M z cos 9 - L x sin0), (11) 

where the x axis is chosen along the direction of the stag- 
gered order. In the self-energy E unp of the unpaired spin 
subsystem we include, for each unpaired spin, the inter- 
action with nearest unpaired neighbors only; this gives 

E unp = 2c(-^Z l Js 2 sin 2 9 - Hscos9). (12) 

Here Z\ is some number which has the sense of an average 
number of "nearest neighbors" for a fictitious random 
lattice formed by unpaired spins, and the energy is again 
taken per dimer, hence the factor 2c in front. The number 
Zi can be safely put to 1, as the following estimate shows: 
We would like to know how many unpaired spins in the 
closest vicinity of a given one (which we put at the origin) 
can have "almost the same" effective interaction J c g with 
it. Since J e g decays as e~ r ^, we have to count impurities 
in a spheric layer with some radius R (the distance to the 
closest impurity, on average R ~ c -1 / 3 ) and thickness 
h r*> £. Thus the average number of other unpaired spins 
which are at the same distance from the origin as the 
closest one is 4nR 2 £c <C 1, so Zi ~ 1 + AircR 2 ^ ~ 1. 

We now derive from the mean-field approach the static 
properties for fixed effective Jnteraction J and afterwards 
perform an averaging over J. It is convenient to param- 
eterize the vectors A and B as follows: 

A = {sin a cos 7, 0, 0}, B = {0, sin a sin 7, 0}. 

It is easy to see that a = 0, 9 = is always an extremum 
of E describing a state with a fully polarized unpaired 
spin subsystem and a non-magnetic (singlet) dimer sub- 
system. A trivial analysis yields the following equation 
for the stability boundary of this solution: 

(H 2 - J 2 )(H - Js) + cs(ZJ') 2 {J + csH - cs 2 J) 

+ZJ'(H - Js)(J -2csH) = 0. (13) 

Another obvious solution is the saturated state with all 
spins fully polarized, which corresponds to 9 = 0, a = 
7r/2, 7 = 7r/4. This state becomes unstable if the field 
drops below the saturation value 

H s = J + ZJ'(l + cs). (14) 

At fixed J' the equation (|13|) has two solutions for H. 
Up to the second order in J', J and to the first order in 
c the lower and upper critical fields are given by 

Hoi ~ Js + cs(ZJ') 2 /J, (15) 
H c2 ~ J - ZJ'(X/2 — cs) — (Z J') 2 /(AJ). (16) 




FIG. 2: Schematic behavior of the staggered magnetization 
as function of the field H. Curves 1, 2, and 3 correspond to 
the regimes a<ii, x\ <Ca<i2, and a S> £2, respectively, 
a = ZJ'/(2J) and jci.a as in JUJ. 

It is easy to obtain the behavior of a, 9, 7 in the vicinity 
of the lower critical field H c i : 

9~ [2(1- H/H cl )} 1/2 , a ~ csZ J' 9 cos 7, 

7 ~ (J — cZJ')s/2. (17) 

The total staggered magnetization per dimer in the vicin- 
ity of H c i vanishes as a square root: 

m st ~ cs(l + ZJ'/J)[8(1 - H/Hci)] 1 ' 2 , (18) 

where the term proportional to J' describes the contri- 
bution of intact dimers. Thus, the following physical pic- 
ture emerges: At zero external field, the presence of im- 
purities induces a finite staggered order, a / 0. When 
the external field H is switched on, it induces canting 
of spins in the field direction, i.e. 9 starts to deviate 
from 7r/2. However, with increasing H the canting angle 
9 moves towards 0, diminishing the effective staggered 
field from the unpaired spins acting on the dimer subsys- 
tem. Unpaired spins become saturated at H — H c \ (9 
becomes zero), and simultaneously the staggered order 
in the dimer subsystem vanishes (a becomes zero). The 
system enters a disordered phase where the dimers are in 
a singlet state and unpaired spins are fully magnetized. 
The resulting phase diagram is sketched in Fig. ^ 

As a byproduct of this description, one can study how 
the high-energy spectrum (above the gap A) is affected 
by presence of impurities. One may replace -Edim in the 
Lagrangian by the full energy E and assume that 
9 is frozen at its equilibrium value. This amounts to 
neglecting the dynamics of the unpaired spiri^ subsystem 
which occurs at a much lower energy scale J, justifying 
to take into account only the static part described by 
a fixed 9. In that way, one obtains that the triplet gap 
splits due to the appearance of the staggered order in the 
x direction; at H — 0, the gaps are related to the order 
parameter A oc m st as follows: 

A y , z = (A 2 +4v 2 Al) 1 / 2 , A x = (A 2 + 12v 2 A 2 ) 1 / 2 . (19) 
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For Ao <C 1, the square of the order parameter deter- 
mining the intensity of the Bragg peak is approximately 
linearly related to the gap energy. This proportionality 
between the change in the gap energy and the intensity of 
the Bragg peak was observed experimentally in TICUCI3 
doped with Mg 01 by varying the temperature. Our 
calculation is for T = 0, but we believe that this relation 
continues to be valid for finite temperature as well. 

Averaging over the impurity distribution. - As final step 
we have to average the staggered magnetization H18(l with 
H c i given by l|15|) . using the distribution function 
Formally, the field Hq at which m st completely disap- 
pears is determined by the upper cutoff of P(J). How- 
ever, the singular (^-function) part of P{ J) has to be 
neglected for two reasons: (a) it has a small weight of 
Zc, so that the corresponding contribution to 77i st is of 
the order of c 2 ; (b) it describes unpaired spins located at 
neighboring sites, and they will have a strong tendency 
to build a singlet state since their interaction with each 
other is much stronger than with the rest of the unpaired 
spin "network" . The regular part P reg of the distribution 
carries the main weight 1 — Zc and corresponds to much 
smaller J with the cutoff at J = Jie~ x ^. Its contribu- 
tion to m s t vanishes above the characteristic field 

i?i = sJie- 1/( - + cs{ZJ') 2 /J (20) 

One has to distinguish two different regimes, depending 
on the strength of the interdimer interaction. For a very 
small interdimer coupling 

^«2a; 1 , where x^e 1 ^ = (2tt) 3 / 2 ^, (21) 

the second term in (|2U|I is always leading, and 

H x ~ H la = cs(ZJ') 2 /J, 

with m s t behaving as a pure square root \fH\ a — H . For 
stronger J' the inequality (|21|l is violated, then the first 
term in l|20|l becomes leading, and 

H x ~ H u = sJie" 1 ^. 

The behavior of m s t for H < Hi depends again on the 
interdimer interaction strength. If (J) re g <C c(ZJ') 2 /J, 
which corresponds to weak interdimer coupling ZJ'/J <C 
2x2 with X2 S> X\ given by 

4 /4 expK 1/2 } = 2(2 7 r) 5 / 2 / C , (22) 

then the main contribution to the staggered magnetiza- 
tion m s t has an overall \JH\ a — H behavior, and addi- 
tionally there is a weak shoulder extending to a higher 
field H = Hu,. The shoulder strength grows with in- 
creasing J' and for 2x2 <C ZJ'/J -C 1 merges with the 
main part, as sketched in Fig. 

In summary, we have shown that in a system of weakly 
coupled spin dimers doped with non-magnetic impurities, 



the application of an external field H may lead to a reen- 
trant behavior of long-range antiferromagnetic order as 
a function of field. From the phase diagram of Fig. 2] it 
is clear that this phenomenon is most likely to occur in 
materials with a sufficiently small ratio of inter- to in- 
tradimer exchange. Therefore, among known materials, 
KCUCI3 is the most promising candidate. 
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